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Abstract 

For a system of N bosons in a 2d liarmonic trap with frequency uj, interacting 

via repulsive forces V <^ Tioj, we develop supersymmetric method to find the 

lowest energy states of rotating Bose condensate as function of two quantum 

numbers, the total angular momentum and the angular momentum of internal 

excitations (generalized yrast states). The energies of these condensed vortex 

states are expressed through the single two-body matrix element of interaction 

V . A broad universality class of the repulsive interactions for which these 

results hold is described by a simple integral condition on V . It includes 

Gaussian, (5-function and /og-Coulomb forces. 

PACS numbers: 03.75.Fi, 32.80.Pj, 67.40.Db, 03.65.Fd 
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Creation of traps for cold rarified atoms |l|], which has made possible experimental in- 
vestigation of Bose-Einstein condensation, has stimulated theoretical studies on the systems 
of weakly interacting bosons confined by parabolic potential Among the important 

questions about behavior of such systems is their response to rotation and the onset of vor- 
ticity In this context, of great importance is the structure and the spectrum of the 

ground states of rotating condensate at given angular momentum, the yrast states [0-^. 

The particular problem of wide recent interest f^-^ arises in the weak coupling limit, 
which is hoped to be reached experimentally in nearest future. One considers 

spinless bosons in spherical 2d harmonic trap 0. The noninteracting system has equidistant 
spectrum hum of high degeneracy which grows exponentially with n at n — ^ 1. This 
degeneracy is related to the number of ways to distribute the total energy hun among the 
particles. The short-range interactions V{r) between the atoms are assumed weak in the 
sense that hoppings between different hun levels can be neglected, and the problem still 
requires nonperturbative solution for the highly degenerate states at single level hun, which 
is similar in spirit to the problem of Landau level for the electrons in magnetic field or to the 
problem of compound states in an atomic nucleus. As is usual for interacting many-body 
problems, evaluation of exact ground state is expensive task. So far, the only results for the 
yrast states in repulsive case were obtained numerically and for simplified 5-forces f^. 

In this work, we provide rigorous analytical solution for the ground state as a function of 
two quantum numbers, the total angular momentum, L{<N), and the angular momentum 
of internal excitations, for a broad class of repulsive interactions. These generalized yrast 
states include the usual yrast states. The results allow transparent physical interpretation. 
In fact, the form of the yrast wave functions which was drawn from numerics for the S- 
interaction case 0, turns out to be valid in general. The universality class of interactions 
is described by an explicit sufficiency condition. 

For this "single Landau level" problem, the Hamiltonian [h = m = 1) is the sum 

H = uHo + V, (1) 



2 



where the first term (the energy of degenerate level) comes from the noninteracting Hamilto- 
nian in the parabolic trap, ^ y + ^rf ^ , and the second (nontrivial) term V" is a properly 
projected interaction V = J2 ^('^y)) see below. In 2d, one uses complex z{z*) = x ±iy in- 

i>j 

stead of vector r = (x, y), and we set = 1 for brevity. With d± = =F ^'§^)^ '^^ employ 
the tetrad of ladder operators a+, a, b'^ and b for each particle, 

a+ = z/2-d+, b+ = z*/2-d^, 
a = (a+)t, b = ib+)\ [a, a+] = [b, b+] = 1, (2) 

which raise (lower) the powers of Zi and z* in the preexponentials of many-body wave 
functions, which are all polynomials times the Gaussian factor |0) = exp{—l/2J2\zk\'^)- 
So, Zi\0) = a,^|0), z*\0) = bj'\0), aj|0) = 0, 6j|0) = 0, etc. The 2d angular momentum is 
the difference L = — L_ of numbers of "up" and "down" quanta, = Y^i ot^k and 
L_ = Yl,i K^k, while the energy of the level in (|l]) is Hq = + + N. The yrast states 
must have minimum energy at given L. Thus, they belong to the subspace with L_ = 0, 
= L, spanned by the homogeneous (degree L) symmetric polynomials , polyg{af)\0) 
which do not involve z*. The basis is formed by polynomials [li,l2, ■■■Jn], obtained by 
applying the operator of symmetrization over N variables, Ps, to the monomials m § 

[h,l2,...,lN]=Psm, m = z['zl\.J^. (3) 

with a given partition of the integer L = J2^n, IK- In our case, L<N, the dimensionality of 

n 

the basis is p{L), "the number of unrestricted partitions" of integer L 0. 

The system enjoys an additional conserved quantity 0,^. The collective operators A'^ = 
Y.f=i and A = {A'^y (with = 1) commute with H, and so does the number of 

collective quanta, The mutual eigenfunctions of H, L and A'^A = v are factorized 

^iL,v) = Zyolyt^Cz,)e-^'-^, (4) 

with Zi=Zi—Z and Z=Y,i Zk/ N. In (^), the total angular momentum L is redistributed be- 
tween internal excitations {J=L—v) and the collective motion (v). It is therefore interesting 
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to consider the ground state as a function of both L and v. We call the eigenvalue of A^A 
seniority for brevity. At given L, the allowed values are w = 0, 1, 2, L — 2, L; the value 
f = L — 1 is forbidden because poly\{zi — Z) =0. The states with definite seniority, of type 
(|), can be obtained by applying the projector Vv to (|^). 

In a given sector L, the projected interaction V in (|I]) is given by pT| ,p 



= PL_=0^'L+=Lm+=LPL_=0. (5) 

In order to evaluate we use the Fourier representation 

oo 



■oo 



where = /q°° rrfr Jo(r|g|)V"(r)/(27r) is the 2d Fourier transform of the potential V{r), 
with Jo the Bessel function [|] and q±={qxr^n.iqy) / \/2. The two-particle combinations 
alj=^{al —Oj) , bfj=^{bf —bj') etc. came from resolving Xi^) and from (H). Sub- 
stituting in and using Baker-Hausdorf relation e"^~" = e' 2' 'e"^e~°, we obtain 



TV 

expansion of V in terms of the two-particle operators B'^ = J2 



]>i 



2[L/2] [L/2] 
fc=0 fc=0 

\4 = - B'^-') + (52fc+2 - s,k+r)B^'^\ (7) 

The highest possible order of B^ is L for L even, and L — 1 for L odd; B^=N{N —1) /2. The 



shape of potential \^(r) is described by the integrals with Kummer function M 0, [13 

<-=^ dt 



Sk= -M{k + 1,1, -t)ViV2i). (8) 
Jo k\ 



Expansion is exact for any interaction V{r) whose moments Sk are finite ||TT],|T2 



Regular methods to obtain exact ground state without solving the whole spectrum are 



not available. We use the approach []I2[ which we loosely nicknamed "supersymmetry" |]14 
Suppose that the Hamiltonian can be written as a sum, 

V = Vo + Vs, (9) 
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such that 

(a) the first term, Vq, is simple, and one can find out its ground state |0) with eigenvalue Smin 
(possibly degenerate). If the second term, Vs, has the two "supersymmetric" properties: 

(b) Vs annihilates the state |0), so Vs\0) = 0, 

(c) Vs is non-negative definite, Vs > 0, (it does not have negative eigenvalues), 

then the state |0) will still be the ground state for the full Hamiltonian Vq + Vs, with the 
same eigenvalue Smin- Indeed, (b) implies that |0) is an eigenstate for the sum Vq + Vs 
with its eigenvalue intact. As one knows form linear algebra, if an Hermitean operator Vq 
is perturbed by a non-negative definite Hermitean operator Vs, the eigenvalues can only 
increase (see, e.g. |1^). This means that the states other than |0) can gain energy [Cf. (b)]. 



As |0) is already the ground state for Vq, it will be the same for Vq + p6| , p!7[ . 

The same arguments apply to the case with additional conserved quantum number v, 
such as [f ,V^o] = b5^5]=0. In the above scheme, the single state |0) is replaced by the set of 
states \L,v), each having the property (a) and (b) in their f -sectors. This is illustracted in 
Fig.l, where the spectrum of Vq is taken degenerate in each f -sector. 

In our case, the supersymmetric triad Vq, Vs and |0) can be established by inspecting 
action of terms Vk in (|^ on the states of partition basis (^). While the first term Vq, by 

virtue of identity NJ^f cif^^i = | Y.fjiO't ~ '^/)(ctj — o-j) + o,t)iJ2i^ cij)-, is reduced to a 
combinations of quantum numbers, 

l/o = {N/2)[{N - l)so -{L- A^A){s^ - s^)], (10) 

the simplest state \L) in (^) with partition [1, 1, ...1] is annihilated by the remainder of the 
Hamiltonian, 

\L) = Psatat...at\0), {V - Vo)\L) = 0, (11) 

Indeed, [0^2, a^ct2^] = -2(ai2ai2 + 1) and thence ai2^^^a]^a^a;]"...aj|0) = and i?''>^|L) = 0. 
From the same commutator, we have [ai2al2,0't'^t]\^) = ^1^02^] |0) = — cti2^|0), and 

thus (5^ — B^)\L) = 0. Substitution Zi = Zi -\- Z [see @)] transforms \L) to a sum 



1^) 



Ps~Zl~Z2-~ZL + ... + Z^-^Ps~Zi~Z2 + Z 



|o) = E^^I^) (12) 



of exactly L seniority observing states of the form each being the eigenvector of A^A 
and therefore of Vq, with the eigenvalue (p^O]). We notice that each f -sector is represented 
by a single term in (^^, identified with Vy\L). 

The spectrum of Vq (see Fig.l.) consists of L equidistant (except 5f(i;)-fold 
degenerate levels with energies ([lOl); g{v)=p{L—v)—p{L—v—l) for v<L—2 and g{L)=l. 



Each f-level contains one and only one state \L,v) = Vv\L) from the sum ([T^ Therefore, 



the set of states \L,v) obey the criterion (a) with the operator Vq. The property (b) with 



V s=V—Vq holds by virtue of (|TT]). In particular, (a) together with (b) mean that |L, v) are 
the eigenvectors of V=Vq+V s with eigenvalues (p!OD. This holds for any interaction V{r). 

The supersymmetric representation for the triad Vq, Vs and \L,v) = Vv\L) would be 
complete if we succeeded to prove non-negative definiteness Vs>0 of the remainder of the 
Hamiltonian V — Vq, criterion (c). So far, we did not specify form of the interaction V{r) in 
(0). We will study now general case and specify the class of potentials which have Vs > 0. 

We have to check signs of all the eigenvalues of Vs in the partition space (^). To avoid 
solving the whole spectrum, in the space of symmetrized states, we use the following trick. 
By definition, the nonzero eigenvalues of Vs in the space (^ coinside with the nonzero 
eigenvalues of PsVsPs in the full space of monomials m in (|^). This latter space has 
dimensionality much higher than p{L) and it includes wave functions of all possible sym- 
metries, including boson sector (fully symmetric), fermion sector (fully antisymmetric) etc. 
In this extended space, the analysis of signs of eigenvalues is however crucially simplified, 
while the contributions from the symmetric sector can be accurately separated. Using the 
symmetry of Vs under permutations of particles, we write 

PsVsPs = PsT.Vs,jPs = iN/2){N-l)PsVs,i2Ps (13) 

where Vs^ij is the contribution from pair of particles i,j to Vs [Cf.(|^,^,^]. In order to see 
that PsVsPs>0, it is sujficient to show that V5,i2>0, because application of any projector Ps 
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from both sides in ([T3|) can add new zero eigenvalues, but can not add negative eigenvalues. 



This follows from the known "inertia theorem" of linear algebra |T5|. Now, we study the 
eigenvalues of Vs^u- Let be the operator of permutation of variables 1 and 2. The triad 

forms a set of mutually commuting operators. Indeed, Vs^u is expressed in terms of B^2 = 
'^12^ '^12 (0)]- [c^i2?'^i2]=l) s-^y can be expressed in terms of 0+12012, using boson cal- 
culus formula ai2'^aj2=Q'i2'^i2(0'i2fli2~l)--- (0,12(112— k+1). The triad T is diagonalized simulta- 
neously right in the basis of monomials m = {z'l Z2 z'^ . . . z'']y } , with the only substitutions 
zi^ ■^{zi—Z2), Z2—^^{zi+Z2). In this basis, the eigenvalues of the triad T depend only on 



1 1 in the subf actor 



■^(2:1— 22) of monomial m, they are {(— h, A^^}, respectively. 



Using Eqs.(|3) and (H) and the summation formula J2k=o k\{N-ky. -^(^ + ~ m ' 



we 



obtain \i, = J rdrV{r)fi, (r) with fi, (r) = e^'^'/^ [#1771 - 1 + T " t)] • The eigenvalue 
(— 1)'^ of 7ri2 helps now to separate out states with wrong symmetry: the eigenvectors with 
li odd are antisymmetric in zi, Z2^ and the projector Ps in ([l3|) kills them all. All even values 
of l\ {<L) can contribute to the bosonic sector, and the corresponding A's must be checked. 

Now, inequality A2„>0 for any 2n<L is the sufficient condition which defines the class of 
potentials for which Vs is non-negative definite, and the triad Vq, Vs = V — Vq and {Vv\L)} 
obeys the supersymmetric representation with properties (a), (b) and (c), with Vy\L) being 
the ground state in its sector L,v with the energy Smin equal to eigenvalue of Vq. 

With this criterion and (^) and ([10|) at hand, we can formulate very general result: For 
any bona fide two-body potential V{r) which satisfies the integral condition 

for any n < L' /2, (14) 



> 



the eigenstates of the Hamiltonian H with minimal energies at given pair of v and 

L(< min{L', A^}) have universal form 

/ o \ N-L+v AT 



with the energies given by the simple moments of V{r), 

S^UL, v) = L + N+ ^^^~^\ o + ^^^^^N{v - L) 

roo ftx) f \ +2\ 

Vo = y^ dte-'V{V2t), Vi = y^ dte-'{- + J\V{^^2i) (16) 

which are equal to expectation values of interaction between two bosons both in the ground 
Vo=(00|l^|00) and the first excited state of oscillator Vi=(ll|l^|ll), respectively pf. 

At fixed L, we have exactly L such equidistant generalized yrast states, marked hy v = 
0, 1, 2, L (t; 7^ L — 1), see Fig.l. Each such state is the "ground state" in the sector L, v 
(of course, there are other states in each sector with higher energies). 

The usual yrast states minimize Smin{L,v) with respect to seniority, v. From ([T6|), they 
can have either f = or f = L. As can be seen, in the domain of validity the first option 
is usually realized. Inequality (Vq — Vi)(f — L) < means that internal rotational excitations 
are energetically favorable, once the interaction energy between two bosons in the state z\0) 
is less, then in the state |0). Physically, the yrast wave functions (^) with v = correspond 
to condensation to a vortex, rotating around the "center-of-mass" , as discussed in The 
maximum seniority states v = L, which correspond to purely collective rotation 0] with no 
internal excitations, were shown P] to be energetically favorable for attractive 5-forces. 

Consider repulsive {Uo>0) Gaussian interaction, V{r) = Uq ^ ; whose range R can be 
varied from zero to infinity. From (p!^, we have A2„ = ,r(2+V) [( 2+^^ )^"— l+4n ^^^^2^ ] > 0, so 
(|1|) holds for any R, and the spectrum is £^in = L + N + ^^^^[N{N -I) /2- ^^^N{L- 
v)]. The yrast states have f = 0. In particular, in the zero range limit, -R^O, we have the 
5-function repulsive interaction V=Uq 5{r). We have A^^ = -^[^^.o + ^ ^ !]• It is seen that 
while Ai=— 1/2, for any li=2n even the condition A2n > (|T^ holds. The energies (0) are 
Ejnin = L + N + ^N(2N — L + V — 2). Yrast states have v = 0, and EminiL, 0) agrees with 
that obtained numerically see also [|T9f . 

For the 2d Coulomb interaction V = Uolog{l/r) with Uq > 0, we have = -^[3?^ — 
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1/n + — 2ilj{2n)] > 0, where ip is digamma function, so (0) holds, and Smin = 

L + N - ^[{log{2) - -f){N - 1) - (3/4)(t; - L)] with 7 = -i^il) = 0.57721... the Euler 
constant. The yrast states correspond to f = 0. 

Let us see, that condition (|1^ imposes only weak restrictions on the repulsive forces 
V{r). Indeed, at small r ~ 0, the factor-function /2n(r) in J^rdrf2n{'r)U{r) in (0) approach 
positive constant values, n — 1 and dfn/dr\r=o = 0, while the first node of the functions f2nij') 
occurs at tq = \j 2\j 12 — 2v^30 ~ 1.43 which is of order of the oscillator length, in our units 
{h=m=uj=l). Therefore, for any short range (as compared to the characteristic length of 
the trap) interaction, the condition ([l^) reduces to 

j d^rV{r) > 0, 

Thus, the results ( p!5| , p!6D hold for short-range interactions, which are repulsive on average. 

The condition ( p!^ holds even for many long-range interactions. This is seen from be- 
havior of function /4(r) (ra = and n = 1 give / = 0) which is positive at r < tq and 
r > Ti = \j 2\l 12 + 2-\/30 ~ 3.10, and is negative only in the interval tq < r < ri {fn{^) 
for higher n behave similarly). It is clear that (p!^ ) holds, if V{r) decreases monotonically 
and fast enough, as is the case for the long range Gaussian and Coulomb forces. One can 
therefore summarize that condition holds for any physically meaningful repulsive interaction. 

To conclude, in the problem of weakly interacting Bose condensate in 2d trap, we devel- 
oped an exact expansion for the interaction to powers of ladder operators. A supersymmetric 
representation for this interaction has been developed to obtain the states with minimum 
energy at given angular momentum and seniority (generalized yrast states, including usual 
yrast states). The energies of these condensed vortex states are given by simple integral mo- 
ments of the interaction potential. A broad universality class of the repulsive interactions 
for which these results are valid is described explicitly. 

Further applications of these results are straightforward. In view of 0, the three- 
dimensional case can be done using the same method. It is also interesting to study region 
of higher angular momenta L > N, where the structure of the basis polynomials will be 
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changed [§], while the numerical studies indicate signs of phase transition The method 
of "supersymmetric decomposition" developed here is not restricted to this particular prob- 
lem and can be applied to fermions and even to the particles with parastatistics. The work 
was supported by FAPESP. 
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Figure Captions 

Fig.l. lUustraction of supersymmetric decomposition. The spectrum of Vq (left) is 
sequence of degenerate levels, labeled by the conserved quantum number v. The spectrum 
oi Vq + Vs is shown on the right. Supersymmetric perturbation Vs>0 splits each level, 
pushing the states up and leaving the lowest energy in each i>-sector intact. 
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